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= m c r y/ I x + ^c r x /*y 
- m c r y /l x 

= n 'c r x/^y 

= moment of inertias in the x, y, z system of 
the space station treated as a rigid body 
= moment of inertias in the x, y, z system 
for body 1 , 

= principal moment of inertias in the x 2 , y 2 , 
Z 2 system for body 2 

= initial product of inertia disturbances for 
the space station. 

= feedback gains 

= torsional stiffness coefficient for cable 
connections' , 

= control mass magnitude 
= undamped crew-quarter wobble frequency 
= undamped torsional natural frequency 
“ (l£ + I 2 )k t /l£l 2 

= Iy/Zy ■ 

= Xy/Xy 
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- coord inafes of control mass m c in 
X/ y, z. coordinate system 

= control force 

a W 

= displacement of the control mass m c 
(see figure 2) . 

■ V 0)/(I * ' V 

■ - V 

= damping factor for the ith mode due to 
feedback control 

/ 2 

- imaginary part of the ith feed-back pole 

- transverse angular velocity components 
of the crew quarters 

= nominal space station spin velocity 

• » 

= e 

= torsionax rotation of body 2 relative to 
body 1 . 
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chapter i 

INTRODUCTION 

1.1 ADMINISTRATIVE INFORMATION 



This report covers the work accomplished on NASA contract 
NASB-27952 during the time period 1 September 1971 to 1 Novem- 
ber 1972 for the George 0, Marshall Space Plight Center, Marshall 
Space Plight Center, Alabama, 35812. The principal investigator 
for this study was Dr. Dara W. Childs, Associate Professor of 
Mechanical Engineering, The University of Louisville, Louisville, 
Kentucky. Inquiries about the results of this study should be 
directed either to Dr. Childs or Mr. Harry Buchanan (S&E-AERO-R) 
Marshall Bpace Plight Center. 


1 . 2 BACKGROUND INFORMATION 

The Principal investigator for this study participated in 
the NASA-ASEB £ ^ mmer faculty program at Marshall Space Plight 
Center during the summer of 1970. The research topic investi- 
gated during this time period was the analysis of a novel "wobble 
damper" for an artificial-g space station. A space craft is said 
to be "wobbling" when its axis of maximum moment of inertia is 
not aligned witn’ its moment-of -momentum vector. This undesirable 
vehicle motion can be controlled by the use of reaction jets, 
passive dampers, control moment gyros (CMG's) and momentum wheels'. 
Each of these devices has both advantages and disadvantages, and 
Dr. Eugene Worley (S&E-AERO-R) suggested the possibility of an 
active system which would employ a single movable mass to generate 
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control torques. The stimulus for this proposition was the 
observation that astronaut motion within the space station 
(moving masses) would be the principal source of spacecraft 
wobble motion. It therefore seemed reasonable to assume that 
controlled motion of a comparable mass could be used to eliminate 
wobble motion. During the course of the NASA-ASEE summer faculty 
program, the correctness of this supposition was verified by a 
combined analytical-simulation study. The results of this study 
are presented in Reference 1, and demonstrate that a movable mass 
controller (MMC) represents an extremely attractive alternative 
for the wobble damping of spacecraft which basically behaves like 
a rigid body. 

The present study was designed to establish the feasibility 
of an MMC in the attitude stabilization of a cable-connected 
artificial-g configuration which can not reasonably be idealized 
as a rigid body. The dynamic model for cable-connected con- 
figurations employed in this study accounts for the aggregate 
motion of the space station and relative torsional motion between 
the crew quarters and counter weight. The development of the 
model is the subject of the following chapter. 


I •itraffftifij 
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CHAPTER IX 
DYNAMIC MODELS 


2 . 1 INTRODUCTION 

Most of the past analyses of cable-connected rotating space 
station configurations have been concerned with (a) the dynamic 
formulation of the problem, or (b) the stability of the system. 
Typical dynamic formulations have considered the extensional 
motion of the system 2 ' 3 or have considered lateral motion of the 
cables 4 ' 5 with the end bodies assumed to be particles, Stabekis 
and Bainum 6 have also examined the problem with an extensible ca- 
ble but consider the end masses to be rigid bodies. 

Posnansky and Heeschen 7 have reported the results of an 
elaborate lumped-parameter simulation of a cable-connected con- 
figuration which includes both lateral cable dynamics and rigid-body 
■* 

descriptions for the end bodies. In terms of the present work, 

their most significant finding was that the fundamental natural 
» 

frequency of the lateral dynamics (coupled cable and rigid-bod^ 
motion) was on the order of 1 cps. By comparison, the rigid-bc^y 
wobble frequency may be on the order of .01 to .1 cps. However, 
since cables provide, comparatively little torsional rigidity, the 
fundamental torsional frequency may well be below the wobble fre- 
quency. From these observed frequency relationships, the dynamic 
representation derived in this study is based on the following two 
assumptions : 

(a) Lateral and extensional motions are a secondary con- 
sideration in designing the attitude-stabilization 
system, and 


kfiiffliwupm ' in -'■’ff.rriM'ii'q 
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(b) Relative roll motion between the end bodies is of 
central importance in the design of th« attitude- 
stabilization system. 

In the sections which follow, a dyamic formulation is derived 
for the space-station/MMC system based on these assumptions. 


~ " - r -pjip^ l fj W l | W.'ii | 1 . l| lll-! L* 'j .' 1 1 1 . J IjJI H— W I I 


2 . 2 GOVERNING EQUATIONS 

Figure I illustmes the basic kinematic variables required to 
define the problem, The X, Y, 2 axes define an inertial coordinate 
system, The x Jt y t , z x and x 2 , y 2 , ? 2 axes are fixed in rigid bodies 
1 and 2, respectively, and their origins coincide with the mass centers 
of these bodies. The x, y, z axes are parallel, respectively, to the 
Xi> y lf z x axes, and the origin of the x, y, z system coincides with 
the net mass center of the system. The x 2 , y 2 , z 2 axes are principal 
axes for body 2. The angle 0 defines the (only) relative motion between 
the x, y, z (xj, y x , z^) axes and the x 2 , y 2 , z 2 axes. The vector R 
locates the origin of the x, y, z system in the X, Y, Z system. The 
angular velocity of the x, y, z system relative to the X, Y, Z system 
is defined by the vector Q, and the angular velocity of the x 2 , y 2 , z 2 
system relative to the x, y, z system is defined by the vector w. 

Body 1 of Figure 1 corresponds to the crew quarters, while body 2 
corresponds to the counterweight. By spinning the system about a 
transverse axis of maximum moment of inertia (nominally the z axis) an 
artlficlal~g environment is induced in the crew quarters. 

The vectors p l and p 2 denote the position vectors in the Xj, y x , z x 
and the x 2 , y 2 , z g systems, and the vectors a 1 and a 2 locate the origins 
of the x 1 , y x , z and x 2 , y 2 , z 2 systems in the x, y, z systems. Hence, 
by definition 






dM mm 

nwh 
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/ p 1 dm « 0 
mi 


/ p 2 dm " 0 
m2 


mi a 1 + mza 2 « 0 


(X) 


where mi and m 2 are the masses of bodies 1 and 7, respectively. 

The following two derivative operations with respect to time are 
required. 


V » 


dV 


dt 


A 

V 


X,Y,Z 


dV 

dt 


x,y,z 


where V is an arbitrary vector. The following notation will be used to 
identify the coordinate system in which the component description of a 
vector Is stated 1 (V) implies a component description in the x, y, z 

or xi, yi , zi system, while (V)^ implies a component description in 
the X 2 , ya, z?. system. These vectors are related by the coordinate 
transformation 


«v ) i2 - [e](v> 


c0 0 -s0 

0 10 
80 0 c0 


(V) 


( 2 ) 


where s0 = sin0, and c0 = cos0 . 

The first equation employed to define the vectors^, to is 

T r = J (a 1 + p 1 ) x (ft + a 1 + p 1 ) dm 


mi 


( 3 ) 


■f J (a 2 + p 2 ) x (R + a 2 + p 2 ) dm 
■ m2 


where T r is the resultant external torque applied to the system at the 

% 

origin of the x, y, z system. By noting the kinematic result 


a* = fi x, a^ + fi x (fix a*) ; i = 1 , 2 


( 4 ) 












and by Introducing Eq. (1), Eq. (3) reduces to 


T r * 2 a* x (ft x m.a i ) + ft x 2 a* x (0 x m,a^) (5) 

1-1 1*1 

+ / (p 1 x p 1 ) dm + / (p 2 x p' 2 ) dm 
mi m2 


By substituting from the kinematic equations 

p 1 «(Qx p 1 ) + ft x (ft x p 1 ) (6) 

A 

p 2 = (ft X p 2 ) + [(w) + (fi x w)] X p 2 
+ (ft + 0)) X [(ft + 0)) X p 2 ], 


implementing the vector identity Ax [B x (S x A)] = B x [A x (B x A)] f 
performing the indicated Integrations, and gathering terns, Eq. (5) 
becomes 


<T r ) * [j](ft) + [(ft)]rj](ft) 

A 

+ [J 2 ](w) + [(w)][J 2 ](w) 

+ [J 2 ](ft X w) + [(w)][J 2 ](ft) + [(ft)][J 2 ](w) 
where the notation [( )] implies 


[(V)] = 



(7) 


and performs the matrix equivalent of the vector cross-product opera- 
tion. Returning to Eq . (7), [j] is the inertia matrix for the entire 
system in the x, y, z coordinate system. Specifically, 


9 


[j] « [J 1 ] + [J 2 ] - l m [(a i )]r(a 1 )] (8) 

i»l 

where [ J 1 ] and [j 2 ] are the inertia matrices for bodies 1 and 2 , 
(relative to their Individual mass centers) in the x, y» z system. The 
matrix [j 2 ] is further defined by 

[j 2 ] - [efcjyt©] (9> 

where [j^] is the inertia matrix of body 2 in terms of the X 2 , y 2 > Z 2 
system. The last three terms in Eq. (7) can be combined as follows 

[(d)](u>) = [J 2 ](fi x w) + |>)][j 2 ](fi) + r(0)]i;j 2 ](u>) (10) 

where the vector (d) is defined by 

(d) = (j[u] - 2[J 2 ]>(fi) , j = trace [j 2 ]. (11) 

Hence, the final form for Eq. (7) is 

(T r ) = [J ](fl) + [(ft) 3 [j 3 (ft) (12) 

+ [J 2 ](u) + [(01)][J 2 ](OJ) + [(d) 3(u) 

The definition of the vectors ft and tu is completed by the inde- 
pendent equation 

T c = / p 2 x (R + a 2 + p 2 ) dm (13) 

m 2 

c 

where T is the resultant torque applied to body 2 at its mass center. 
Substituting from Eqs . (1) and (6) into Eq. *(13), integrating over m 2 , 
and then substituting from (10) and (11) yields 
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(x c ) i2 - [jyoo t , + C(® l2 ][j| 2 ](fi) 12 w> 

+ + [w l2 KJi t 3<“>i 2 + 


Since the X 2 , yz , zg axes are principal axes for body 2, the y component 
of Eq. (14) yields 

T y = l y + l2 y ® + (I 22 ' - npoe-tf - n x n z c(20)} (15) 

2 2 2 

where 1 1,1 , I are the moments of inertia for body 2 relative to the 

X 2 y 22 

x 2 > Y2» z 2 Principal axes, respectively. Further, ft , ft , ft are the 

x y z 

components of ft in the x, y, z system. Eq. (15) and the matrix equation 
(12) define the uncontrolled motion of the system. The assumption is made 


in this study that is the reaction torque defined by 

T C = - k 0 

y t 

, 1 

where is a torsional spring constant. 


(16) 


Figure 2 illustrates the MMC system to be used ip this study. The 
essential element of the stabilizer is the point mass m^ whose motion 
is limited by a tube (or other physical constraints) to be parallel 
with the z axis. The tube is attached to body 1, and the position of the 
mass within the tube is defined by the variable w. The governing equation 
of motion for w is 3 


• 

f = m {W + (ft r -ftr) + rftft 
z c L x y y x xf z 2 


(17) 


+ r ft ft - w(ft 2 + ft 2 )} 

, y z y x y J 

where is the z component of the reaction force apt Lied- to m^, and 


r | 
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.(r , r , w) are the components of r, the vector which 

x y 

locates the control mass in the x, y, z system. 

The attachment of the mmc system to body 1 as illustrat- 
ed in Fig. 2 will only change Eq. (12) . The required modifi- 
cation is 

(T r ) = |7] (si) + [( Si )] [J] (S2) (18) 

+ (J 2 ] (w) + [(w>) [J 2 ] M + [(d)] M 

A 

+ m c f(r)] (r) + 2m. c [ (r) ] [ (ft) ] (r) . 

where 

[J] = [J] - m c [ (r) ] ( (r) ] (19) 


Summarizing, the variables (ft , ft , ft_,0, w) are defined by 

x y z 

Eqs . (15), (17.), and (18). 

From Eq. (19) , the inertia matrix [J) is seen to be a 
function of both w and e * Of particular importance are the 
relationships 


Z xz = W 0 ’ + m o r x w - (I Z \ - X xz )s8ce 


J yz m V (0) + “o r y w 


2.3 UNCONTROLLED MOTION: LINEARIZED MODEL 


The free motion dynamics problem can be summarized as 

follows. Undei 'ideal conditions, the H (moment-of-momentum) 

vector and the axis of maximum moment of inertia of the 

\ 

spacestation coincide. If (due to either external torques 
or internal mass motion) the two are forced apart, a 
condition generally called "wobbling" results. This motion 
is characterized by a nominally constant rotation rate about 
the s axis accompanied by oscillations in and 9. 
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The defining equations for free motion are Eqs. (12) and 
(15), and a valid linearization of these equations is 
possible based on the following observations: 

(a) ft « ft, the average spin velocity, 

(b) ft and ft are very small in comparison to ft, and 

x y 

(c) 0 is generally small justifying the approxima- 
tions c0 = 1, s0 » 0. 

The resultant linearized equations are 

ft x - - ft^aft + 0e(l~b 2 )ft - aft 2 a(0) 

n y = R x bin + (k t /x^ ) e - bs? 2 e (o)/q, (21) 

0 = - ft (b-ba) ft/qi - (P?+b 2 ft 2 )6 + bft 2 3(0)/qi 

A 


where the result, bi = (b-q 2 b 2 )/qi, is employed . Two 
special cases can be obtained from these general results. 
First, if -»■ “,-0 and 0 ->• 0, qi -> 1, and bi •> b. The 
result is the rigid body equations. 


ft =-ftaft-aft 2 a(0) 
x y ' 

oiy = ft^bft - bft 2 £(0) 


( 22 ) 


with the rigid-body wobble frequency 

\/ n __ 

p == ( ab ) ' * ft (23) 

The second special case results if k,=b 2 =0. Since e - I 2 /I 
■ t y x 

is in general quite small, the resultant .motion is approxi- 
mately defined by 



T 


T 
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b -fl aE - (Ah7 2 c* (0) 

- bE*0 (0//q A (24) 

0 « -jybjE - bJT 2 © + (0)/q 

The "rigid body" wobble frequency and the torsional oscillation 
frequency p 2 are defined by 

„ 1 / 2 _ 1/2 

p : = {a b x } to , p 2 = b 2 (25) 

2.4 CONTROLLED MOTION : LINEARIZED MODEL 

As stated previously, motion of the complete system is de- 
fined by Eqs. t (15),- (17),- and (18). A valid linearization of the 

controlled system is possible based on the preceding assumptions 

* 

for free motion (fi = <<s7;s0 = 0,c0 = l) plus the 

z x y 

following assumptions for the MMC : 

(a) m w 2 /I ccl, m w 2 /I «1, and 

w a y 

(b) the "coriolis-torques" 2m wwfi and 2m ww ft 

c x c y 

which arise in Eq. (18) due to the term 
2m c ( (r) ] [ (ft) ] (r) are negligible. 

The resultant linearized equations may be written as 
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V*y 


-^(b-b 2 )/q l 0 


0 0 


jTr x (b x -l) ftr (a-1) 4Tr e(l-b 2 ) r x k fc /Xy 0 ft 2 e x v 


1 0 



-aa (0) 


-b3 (0)/q x 

> u + ft 2 

b3(0)/q x 


0 

c 

r y aa(0) -r x b3 (0)/q x 


0 


•r /I 
y x 

r /X 1 
x / y 

•r /X 1 
x' y 


(1+e i )/m 


In the equation above, u is the control contribution to 
the force f appearing in Eq. (17) . A comparison between the 
solutions from this linear model with those of the general non- 
linear model (Eqs . (18), (17), and (15) ) for a wide range of 
cases revealed no significant differences, In most cases no 
discernible differences were evident in plots of the solutions. 
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DATA/ AGREEMENTS , AND ASSUMPTIONS 

3.1 CONTROL OBJECTIVES AND CONTRAXNTS 

The general attitude-stabilization control objectives may 
be defined in terms of the control force u and the state-vari- 
able vector X (q x / 3, 0, v, w). Xn terms of these variables, 
the basic requirement is the synthesis of a feedback control 
logic u = u(X) which satisfies the following descending order of 
objectives : 

(a) The attitude-stabilization system is required to eliminate 
the oscillations in the transverse angular velocity components 
of body 1 (n , o ) , i.e., the desired terminal state is 

X v 

fijjttj;) — ( t^ ) — 0. t- 

(b) While the comfort of the astronauts does not explicitly 
depend on 0 also being forced to zero, the coupling of 

the system is such that condition (a) cannot generally (k t ^0) 
be otherwise satisfied. Hence, an additional desired ter- 
minal state is e(t^) = 0. 

(c) It is desirable that the full capacity of the controller 
be restored after a disturbance has been eliminated, hence, 
the additional desired terminal states w(t^) = w(t^) = 0. 

In summary, the desired terminal state is 

^(t f ) = ^(t f ) - e(t f ) = w(t f ) = w(t f ) = 0 ( 27 ) 

The control constraints can be summarized as follows: 

(a) Motion of the control mass is amplitude limited by physical 
constraints, i.e., 


( 28 ) 


| v? | 2 w “ 16,5 ft 
(b) A control which requires an .excessive control force u is 
deemed to be impractical. Specifically, in this study 
control forces in excess of thirty-five pounds were deemed to 
be unacceptable; hence, 

| u | 5 U « 35 lbs. (29) 

In addition to these specific and easily quantifiable 
requirements and constraints, the feedback control logic to 
be derived was constrained to be time-invariant and (generally) 

linear in form. Specifically, a control law of the form 

• • t 

u - ka S7 + k 2 ^,, + k 3 6 + k t* 6 + k$w + Kgw + k 7 f w(T)dt (30) 

x y 0 

was, employed where the integral term was added for the sole 
purpose of enforcing the terminal boundary condition w(t f ) = 0, end 
has no other intended influence on the dynamics of the system. 

From Eqs ’. (26) and (30), the control system has fche form 


X = AX + bu + £ , u = k T X (31) 

or in closed-loop form 

X = [A + bk T ] X + f (32) 

In general terms, the control synthesis problem addressed here 
can be stated as follows. Determine those gains k^ which will 
cause the controlled system defined by Eq. (32) to rapidly re- 
cover from an initial disturbance, and to approach the desired 
terminal state given in Eq. (27) without violating the constraints 
of Eqs. (28) and (29). 
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3.2 SPACE- STATION CONFIGURATIONS AND DATA 

To this investigator's knowledge, NASA's more recent investi- 
gation of artific.ial-g space stations have concentrated on the 
following three concepts: 

(a) A nuclear-powered, nominally axisymmotric configuration for 
which two bodies are connected by a comparatively rigid con- 
necting tunnel, 

(b) A nuclear-powered, nominally axisymmotric configuration for 
which two bodies are connected by cables, and 

(c) A solar-powered, axymmetric cable-connected configuration. 

Concepts (a) and {b ) were investigated by McDonnell-Douglas Astro- 
nautics company under the direction of Marshall Space Flight Center, 
while concept (c) was investigated by North American-Rockwell under 
the direction of The Manned Space Flight Cancer, NASA Houston. 

* The superior performance of an MMC for uctitude stabilization 
cf concept (a) has previously been demonstrated in Reference 1, 
and the investigations of this study are restricted to concept (b) . 
Specifically, the configuration investigated in this study is 
based on the geometric and inertial properties arrived at by 
McDonnell-Douglas Astronautics in Reference 8. Table A illustrates 
the significant longitudinal dimensions of this space station de- 
sign in its various stages of deployment. 

The inertia properties for the sp^ace station in its various 
stages of deployment are provided in Table B. The asymmetry which j; 

is evident (i.e., l z >I X , 1 l z >.Tx , l| >lj.) in this table was added 

i. 

by the investigators to the nominal symmetric properties provided 
by Reference 8. The following guidelines were used to arrive at 
the^e modified inertial properties: 
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TABLE A; AXIAL SPACE-STATION DIMENSIONS 


body 1 body 2 



Deployment 


Cable 

Length 


0 57.16 166.2 0.0 -51.5 
1/4 69.28 201.4 35.3 -63.6 
1/2 81.41 236.7 70.5 - 75.8 
3/4 93.53 272.0 105.8 - 87.9 

1 105.7 307.2 141.0 -100.0 


Cj (ft) m distance from C.M, of body 1 to C.M. of space 
station 

C 2 (ft) « distance from C.M, of body 1 to C.M. of body 2 
r^(ft) - distance from C.M. of space station to MMC 

(a) Dynamic stability of the system requires asymmetry. 

(b) The crew quarters are -nominally axi symmetric , and the pro- 
vision of asymmetry requires additional design effort and 
expense. The assumption was made that asymmetry would be 
provided by the internal mass-distribution design cf the 
crew quarters, and the asymmetry used in Reference 1 was 
taken to be representative of the degree of asymmetry which 
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T 
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could be obtained by this approach. 

(c) The counterweight is also nominally axisymmetric, and it 
turns out that (from a controls viewpoint) a perfectly 
axisymmetric body is ideal. Consequently, sufficient asym- 
metry was provided in body 2 to account for residual (but 
unintentional) counterweight asymmetry. 

3.3 REFERENCE DISTURBANCES AND CONTROLLER SIZING 

The reference disturbance used throughout this study consists 
of four astronauts (773 lbs. or 24 slugs) moved instantaneously 
to the worst possible locations within the crew quarters. The 
resultant initial product-of-inertia disturbances for the fully 
deployed configuration are 

I (0) = 3250 slug-ft 2 , I (0 ) * 51,000 slug-ft 2 (33) 

x z y z 

* ** 

These product-of -inertias were used for all stages of deployment. 

Eq. (20) demonstrates the linear dependence of the X and 
I product of inertias upon the controller position w, and it 
is primarily through this dependence that the MMC system is used 
to generate control torques. The MMC considered in this study is 
defined by the parameters 

r = 16.5 f , rn ~ 11 slugs (34) 

x c 

The parameter r^ depends upon the stage of deployment and is de- 
fined in Table A. The product-of-inertia capacity of the con- 
troller is from Eqs . (20) and ^28) I C XZ " 3000 slug-ft 2 , X°y Z - 
18,200 slug-ft 2 (r^ = 100 ft) . The equivalent static torque capa- 
bility (for ST = 4 rpm, and r^ = 100 ft.) is approximately l r - c x i, Tnax “ 

I I c IfT 2 = 3180 ft-lbs , and It I = I I c |0 2 = 525 ft-lbs. Since 
1 yz 1 , *y ' max 1 xz 1 








TABLE B: INERTIA PROPERTIES AND PARAMETERS 


Body 1 Body 2 



Deployment 

Stage 

I v (slug-ft 2 ) 

Iv (slug- 

-ft 2 ) 

a 

0 

64.926x 10 6 

65.009x 

10 6 

.9796 

1/4 

90.166x 10 6 

90.249x 

10 6 

.9853 

1/2 

119.91 x 10 6 

119.99 x 

10 6 

.9889 

3/4 

154.13 x 10 6 

154.21 x 

It-* 

O 

m 

.9914 

1 

187.10 x 10 6 

187.18 x 

10 6 

.9929 


Constant Inertial Parameters 


I 1 = . 621x 10 D slug-ft 1 


I 2 = .787* 10 6 slug-ft 2 

y 


I y = 1.407x 10 6 slug-ft 2 


b = .0500 , b l = .0709 , b = .0592 


the products of inertias are linear functions of w, the control 
torques may be varied linearly from zero to full capacity. 


3.4 DEPLOYMENT AND RETRACTION PROCEDURES 

The assumption was made in this study that the nominal design 
spin velocity SI is 4 2:pm, and that in both the deployment and 
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retraction pnases the spin- ‘Velocity could only exceed this value. 

In other words , the rigidly assembled space station would be over 
spun prior to deployment, As deployment proceeded, the spin- 
velocity would drop (conservation of angular momentum) ,and when 
fully deployed would rotate at 4 rpm. Basically the same pro- 
cedure ’would be followed in the retraction phase. Hence, spin 
velocities of 4 rpm and higher were investigated. This approach was 
found to be conservative, since in general increasing the spin- 
velocity degrades the controller's performance. 

3.5 TORSIONAL STIFFNESS COEFFICIENTS 

The two following configurations were used in defining the 

torsional stiffness k . : 

_ .. ... .t _ . 

(a) a zero torsional stiffness configuration (k. t = 0) , corresponding 
to a single cable design, and 

(b) a maximum torsional stiffness configuration corresponding 
to the McDonnell-Douglas design. 

The relatively high torsional stiffness of the McDonnell-Douglas 
design was verified by simulations (at MSFC)of alternative designs. 
The calculation of the values for k^. were based on the work of 
D. Nixon 9 (S&E-AERO-R) , and were aided by personal communications 
with Mr. Nixon. 

3.6 SYNTHESIS PROCEDURE 

One important feature of the state variable formulation 
of a linear feedback control system is the fact that the poles 
of the closed-loop - transfer function can be positioned arbi- 
trarily by the proper choice of feedback gains (if all the state 
/ ' . 

m m «> ; ; ? - 1 1 1 
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variables arc available) , The closed loop characteristic 
equation of a system written in the form X = [A+bk^JX is given 
by the determinant | si - A •« bk“ I . For an nth order system, 
the characteristic equation is an nth order polynomial in the 
Laplacian operator s, with n + 1 terms. The coefficient of the 
s n term equals one, and the coefficients of the s 11 " 1 to s° terms 
provide n algebraic equations in terms of the n unknown feed- 
back gains k^ that can be equated to the corresponding co- 
efficients of any v'.a sired nth order characteristic polynomial. 

By solving this set of n equations, the closed- loop poles of the 
system can be positioned to obtain desired frequency domain 
solutions . 

The effectiveness of the root specification method as a con- 
"tror system design ""technique is limited by the fact that the 
transient magnitudes of the various system and control variables 
cannot be predicted from frequency-domain solutions. Two major 
constraints imposed on the performance of the control system 
considered in this study were a maximum control mass deflection 
magnitude, (Eq. (28)) and a maximum control force magnitude 
(Eq. (29)). While the ability to calculate feedback gains to 
yield desired degrees of damping on the system variables is of 
considerable use, intuitions based on past experiences and trial 
and error procedures provide the only insight into the effect of 
pole locations upon the maximum' variable magnitudes in time domain 
solutions. In this phase of the study, gains were first calculated 
to yield roots specified from physical and intuitive bases. Time 
responses using these gains were then obtained by numerical inte- 


gration of the system's equations of motion to determine 
whether the gains caused the system constraints to be violated. 

The method of closed-loop pole specification was applied 
to both the fourth-order, non-deployed rigid body configuration 
and the sixth-order, cable connected configuration of Eq. (26). 

The characteristic equations of the two systems were established 
by expanding the determinant of the closed-loop system coefficient 
matrix. The closed-loop characteristic polynomial for the cable- 
connected system is provided in Appendix A. Sets of equations 
relating the unknown feedback gains to the desired closed-loop 
poles were obtained by equating the coefficients of the closed- 
loop polynomial to those of the desired characteristic equation. 

A computer program utilizing the IBM Scientific Subroutine Package 
program, SIMQ, was developed to solve the equation sets. Two 
additional IBM subroutines, HSBG and ATEIG, were used to calculate 
the eigenvalues of the closed-loop system matrix for each set 
of calculated gains to verify that the desired roots were actually 
obtained. 

The general criteria for selecting desired root locations can 
be explained as follows. First, one notes from Eq. (26) that 
without control the control mass motion is itself unstable due to 
the coefficient A 56 = Sl 2 e 1 . Assuming that the controller is 
stabilized (via the k c w term in Eq. (30) with all other k. = G), 
the system defined by Eq. (26) has the following three undamped 
modes of motion and natural frequencies : 

(a) wobble motion involving and (natural frequency = P a ) , 

(b) torsional motion involving $ and 0 (natural frequency = P 2 ) , 
and ’ 
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(c) control mass motion involving v and w (natural frequency 
= P 3 ) • 

Hence, the characteristic equation for the uncontrolled system 
has the form 

3 

ir (s 2 + p, 2 ) =0 (35) 

1-1 

The gains are selected so that the closed-loop characteristic 
equation has the form . 

3 

(s + p 0 ) IT (s 2 + 2 CiPfS + Pi 2 ) - o (36) 

i=l 

where the s + p 0 term arises due to the integral feedback term. 

The damping factor for wobble motion was selected to be unity 
in all cases. The damping factors for the control mass and 
torsional motion were selected to yield equivalent damping constants, 
i.e.,? l p 1 = c 2 p 2 = C 3 p 3 < The natural frequency of the control 
mass p 3 and the pole location due to integral feedback p 0 were 
selected in accordance with various criteria depending upon the 
configuration. 
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CHAPTER IV 
RESULTS 


4 . 1 INTRODUCTION 


As explained in the preceding chapter, the synthesis procedure 
employed in this study consists of (a) selecting (guessing) 
closed-loop, root locations, (b) calculating gains k^ which will 
yield these roots, and (c) obtaining the transient solutions to 
determine whether the force and deflection constraints have been 
violated. The trial-and-error character of this procedure is 
such that one is never quite sure that yet another try would not 
improve things. However, the solutions presented here are the 

result of ? large number of such trial and error refinements.,.. ...... 

While some slight improvement might be achieved by additional 
effort, these results in general represent best-possible con- 
troller performa* ze , i.e., optimum selection of gains . Fur- 
thermore, the controller performance through a wide-range of 
configurations and operating conditions is consistently excellent. 


The results presented in the following sections demonstrate the 
influence of the following factors on controller performance: 

(a) stage of deployment, 

(b) nominal spin velocity n, and ' 

(c) torsional , stiffness k t . 

The results consist of both tabular data and illustrations of 
transient solutions. The tabular data (which is presented for 
each case) contains the following information: 
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(a) « ratio of effective to critical damping for a 3 1 modes, 

(b) p » wobble frequency, 

(c) p 2 = torsional natural frequency, 

(d) specified closed-loop poles, 

(e) calculated gains k^, 

(f) peak control force |u| , and 

(g) peak control mass deflection |w| . 

Transient solution illustrations are provided for the two fully- 
deployed configurations (k t = 0, k t / 0) both with and without 
control. In addition, transient solution illustrations are pro- 
vided for the 1/2 deployed configuration (k^ / 0) at maximum spin 
velocity. These illustrations were included to provide a quali- 
tative indication of the controller's effectiveness. 

4.2 UNDEPLOYED SPACE STATION 

The requisite input data for the undeploye^ space station is 
provided in the preceding chapter. As previously noted in Section 
3.4, the nominal spin-velocity can equal or exceed 4 rpm. The two 
cases considered for the non-deployed configuration are (a) 0, = 4 
rpm and (b) H = 11.5 * Dm. An initial value of ST = 11.5 rpm. is 
sufficient to yield a spin-velocity of 4 rpm. for the fully-deployed 
configuration. The results for these two cases are illustrated 
in Tables 4.2(a) and (b) . Inspection of these tables reveals that 
the lower spin-velocity case is the more easily controlled in the 
sense that a higher control force is required. 
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TABLE 4.2(a) ; Undeployed Configuration, 52 = 4 rpm. 

Input Data; 52 = .413 rad./ser, , , pj =-• .100^ rad. /sec. 

Control Data; C] - 1.0, C 3 = .1667 

kj = 1.815 x 10 3 lb. sec. /rad. , k 2 = ‘-3.767 x 10 4 lb. sec. /rad. 

k 5 = .1037 lb. sec, /ft., k 6 = -4.456 lb. /ft. , k 7 = -.1109 lb. /ft. sec. 

Closed-, Loop Pole Location 

f2 x , S2y w f w 

real -.1004 -.1004 -.0288 

imag. 0.0 .5942 

Transient Performance; 

I u I max “ ^.9 lbs, l w lniax = ^ ft. 

‘TABLE ‘4‘. 2 (b) Und'eployed Configuration, 52 = 11.5 rpm. 

Input Data; 52 = 1. 2035 rad. /sec. , p x = . 2891 red./s.ec. 

Control Data: \ x = 1.0, ? 3 = .40 

kj - -1.081 x 10 3 lb. sec. /rad. , k 2 = -1.626 x 10 4 lb. sec. /rad. 

k g = -10.91 lb. sec. /ft., k 6 = -9.965 lb. /ft. , k ? = -.1596 lb. /ft. sec 


Closed-Loop 

Pole 

Location 





Uu f 

w 

/w 


real 

-.2892 

-.2892 

-.0288 


imag 

0.0 

.6626 


Transient Performance 



1 u Imax 

24.4 

lbs . , 1 w 1 

' 1 'max 

= 10.8 

lbs . 
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4.3 CONTROL FOR k t - 0 . 



Since a space station will normally be in the £ ally- deployed state 
the effectiveness of the MMC in providing attitude stabilization 
for this particular configuration is of fundamental interest, 
Evidence for the controller's superior performance is provided by 
a comparison of the uncontrolled motion illustrated in figure 3 
with the controlled motion in figure 4. In addition, from Table 
4.3(a), the peak force and deflection magnitudes are seen to be 
|wl „ = 16.5 ft., and |uL„„ = 14.2 lbs. Note is made that the 
gains on 0 and 0 (k 3 and k 4 ) are zero, and this is characteristic 
of all k^ = 0 cases. Since no disturbance torque is applied to 
the crew quarters due to 0 motion, there is no reason to attempt 
control of the relative torsional motion. 


TABLE 4.3(a). The Fully Deployed Configuration, ft = 4 rpm. , k^ = 0 . 
Input Data: Pj/Sl = , 230 , p 2 /& = .223 
Control Data: £ x = 1.0 , 0., ? 3 = .20 

kj - 610.1 lb. sec. /rad. , k 2 = -1.561 x 10 4 lb . sec ./rad. , k 3 = k 4 = 0., 
k s = -.634 lb. sec. /ft. , k g = -3.974 lb. /ft. , k ? = -. 0910 . lb ./ft. sec. 
Closed-Loop Pole Location 

ft x , fty w 0 / w 

real -.1109 -.1109 0. -.0287 

imag .0093 .5431 .0933 

Transient Performance : 

|u| = 14.2 lb., |w| = 16.5 ft. 

1 1 max 1 1 max 




1 
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4. a Controlled Motion: ft versus t for the fully 

deployed configuration with - 0 , and fl=4rpm 
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4.c Controlled Motion: w versus t for the fully 
deployed configuration with kt= 0, and 
ft = 4rpm. 
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CONTROL FORCE (LBS) 



b . 3/4 Deployed Configuration , ~ ~ 4 and 4.85 rpm 

The results for this case are summarized in Table 4.3(b). A 

comparison of these results with those or Table 4.3 'a) reveals scant 

differences between the 3/4 and the fully deployed configuration 

when H ■« 4 rpm. This comparative insensitivity to the state of 

deployment can be explained by noting that the only factors which 

vary in the governing Eg. (26) due to deployment are a, r y/i x / and 

e = I 2 /I . The latter two parameters are small, and have a com- 
y x 

paratively minor influence on the systems’ behavior. Conversely, 
the dependency of the wobble frequency p^ upon a (Eq. (25)) makes 
it a significant parameter. However, Table B reveals that a only 
changes slightly from the non-deployed to the fully deployed con- 
figuration . 

A comparison of the results of case 1 and 2 in Table 4.3(b) 
demonstrates the marked dependence of controller performance upon 
the nominal spin velocity. The natural frequencies pj , p 2 display 

I 

the predicted (Eq. (25) ) linear dependence upon <“T while the peak 
force l u l max increases markedly with increasing jjT. 

c . 1/2 Deployed Configuration, = 4 and 6.23 rpm. 

The results for these two cases are presented in Table 4.3(c), and 
continue to demonstrate (a) the comparative insensitivity to the 
stage of deployment, and (b) the marked dependence of controller 
performance upon the nominal spin velocity 

d . 1/4 Deployed Configuration, H = 4 and 8, i . pm 

The results for these two cases are presented in Table 4 . 3 (d) , and 
represent basically the same controller behavior cited previously 
for the 1/2 deployed configuration. 
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TABLE 4.3(b): 3/4 Deployed Configuration, k t = 0.0 
Case 1 : il = 4 rpm - .418 rad. /sec. 

Input Data: pj /IT = . 230 , p 2 /IT = .223 
Control Data: t; 1 = 1.0, C 2 ~ 0.,? 3 - .1667 

kj = 1.023 x 10 3 lb. sec. /rad. , k 2 - -2.259 x 10 4 lb. sec. /rad. 
k 3 = k 4 = 0, k 5 = 1.399 lb. sec. /ft., k 6 = -5.3940 lb. /ft. 
k 7 = -1.309 lb. /ft. sec. 

Closed-Loop Pole Location: 



ft x , fty 

V7 

0 

/ W : 

real 

-.1109 

-.1109 

0 . 

-.0286 

imag. 

.0108 

.6554 

.0932 



Transient Performance: 

lul = 19.1 lb., Iwl = 16.5 ft. 

1 1 max ' 1 1 max 

C ase 2 : ft = 4.85 rpm = .507 rad. /sec. 

Input Data: p 1 /~ = .230 , p 2 /ft = .223 
'Control Data: 1.0 , ? 2 = 0.0, £ 3 = .1667 

k x = 914.1 lb. sec. /rad. , k 2 = -2.729 x 10 4 lb . sec. /rad. 
k 3 = k 4 = 0, k 5 = 1.732 lb. sec. /ft., k 6 = -7.947 lb. /ft.. 
k 7 = -.19257 lb. /ft. sec. 

Closed-Loop Pole Location: 


ft x , fty w 0 /w 


real -.1344 

-.1344 

0. 

-.0286 

imag. .0125 

. 7950 

.1131 


Transient Performance: 




= 28. 0 lbs. , Iwl 
1 1 max 'ii max 

= 16.5 ft. 
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TABLE 4.3(c): 1/2 Deployed Configuration, k^ “ 0.0 
Case 1 : J2 - 4 rpm - .418 red. /sec. 

Input Data: p j/fT = . 230 , p 2 /ft = .223 
Control Data: £ x ~ 1.0, C 2 83 0.0 , C 3 = .250 

k L - 268.26 lb. sec. /rad. , k 2 = ”9 . 785 x 10 3 lb, sec. /rad. 
k 3 ~ k 4 - 0, k 5 = -2.3147 lb. sec. /ft., k 6 - “2.794 lb. /ft. 
k 7 ~ -.0580 lb. /ft. sec. 

Closed-Loop Pole Location: 


^x f 


w 


0 


/ w 


real 


imag. 


-.1107 -.1107 0.0 

.0117 .4285 .0932 


Transient Performance: 


-.0286 


u max = 9*7 lb. , w max - 16.1 ft. 


Case 2: Ti = 6.23 rpm = .652 rad. /sec. 

Input Data: p x /fi = .243 , p 2 /fi = .223 
Control Data: 1.0 , ? 2 = 0.0 , ? 3 = .20 

k x = 358.6 lb . sec. /rad. , k 2 - -2.393 x 10 4 lb. sec. /rad. 
k 3 = k 4 = 0. , k 5 = -.9139 lb. sec. /ft. 
k G = -9.619 lb, /ft. , k 7 = -.2206 lb*/ft.sec 
Closed-Loop Pole Locations : 

ft x , fiy w 0 / w 

real -.1726 -.1726 0.0 -.0287 

iniag. .0170 . 8455 . 1453 

Transient Performance: 

i u lmax 5=5 ' ^4.0 lbs., |w| max = 16.5 ft. 


TABLE 4.3(d); 1/4 Deployed Configuration/ k^ ~ 0.0 

Case 1 : ft = 4 rpm = .416 xaci./oec. 

Input Data: Pj/M « . 230 , p ? /P. = .223 

Control Data:?j« 1.0 , £ 2 =s 0.0, .1607 

k x - 1.013 x 10 4 lb.sec./rad. , k 2 = -2.238 x 10 4 lb. sec. /rad, 
k g « k 4 =s 0, k 5 1. 350 lb.sec./ft./ k 6 = 5.362 lb. /ft. 
k 7 as -.1301 lb. /ft. sec. 

Closed-Loop Pole Locations : 



^X t 

w 

0 

/ w 

real 

-.1105 

-.3105 

0.0 

-.0285 

imag. 

.0140 

.65342 

.0931 



Transient Performance 


u = 19.0 lb. , wl „ = 16.5 ft. 

- — max - — — • — ' - **- max 


Case 2 : 8.3, rpm = . 867 rad. /sec. 

Input Data: Pj/fF = . 254 , p 2 /fi = .223 

'Control Data: = 1.0 , £ 2 = 0. , £ 3 = .333 

k : as -357.9 lb. sec. /rad. , k 2 = -1. 060 lb. sec. /rad. 

k 3 = k 4 = 0. , k 5 ••= -7.258 lb.sec./ft., k 6 = -7.880 lb. /ft. 

k ? = -.1400 lb. /ft. sec. 

Closed-Loop Pole Location: 

* 'fly, fly w 0 /w 

real -.2293 -.2289 0.0 -.0287 


-.2293 -.2289 0.0 -.0287 

.0221 .6478 .1931 


imag* 

Transient Performance: 


u ^ = 23.9 lb., w = 12.9 ft. 

max max 
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4.4 CONTROL FOR FINITE TORSIONAL STIFFNESS 
a • The Fully De p loyed Confi q u rn i- ion , JT_ - 4 rpm 

The obvious observation that: a cable-connected space station will 
be fully-deployed throughout most of its lifetime makes establish- 
ment of effective control for this particular stage of fundamental 
importance. Evidence for the controllers’ superior performance 
is provided by a comparison of the uncontrolled motion illustrated 
in figure 5 with the controlled motion illustrated in figure 6. 

From table 4.4(a), the controlled motion of figure 5 corresponds to 
critical wobble damping rates, and results in the peak force 
magnitude |u| max = 16.1 lb., and peak deflection magnitude I w I max = . 

16.3 ft. 

TABLE 4.4(a): The Fully, Deployed, Configuration , ?2 = 4 rpm 

Input Data: k t = 1.463 x 10 5 ft . lb;/rad. 

p x /fi = . 230 , p 2 /fi ss 1.36 

Control Data: Z l = 1 . / K 2 ~ .169, C 3 = .167 

kj = 3061. lb. sec. /rad. , k 2 = -5.063 x 10 ** lb . sec . /rad . , 
k 3 = -3.052 x 10 4 lb * sec. /rad. , k 4 = 2.215 x 10 3 lb./rad. 
k s = -1.764 lb. sec. /ft., k 6 = -5.653 lb. /ft. , k 7 = -.161 lb. /ft. sec. 
Closed-Loop Pole Location: 

, fiy W 0 f W 

real -.1109 -.1109 -.1109 -.0288 

imag . 0 . 0 . 6464 .- . 6559 

Transient Performance : 

Hmax = 16 - 1 lb -' l w lmax “ 16 ' 3 ft ‘ 
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6.d Controlled Motion: w versus t for the fully 
deployed configuration with = 1.463 x 10 s 
ft. lb, /rad., and n = 4 rpm. 



to • 3/4 De p loyed Configurn t ion . U ^ 4 and 4.85 rpm. 

The results tor this case me summarized in Table 4.4(b). A com- 
parison of these results withe those of 4.4(a) for sT « 4 rpm 
reveals that the only significant difference caused by a 1/4 
retraction is an increase in k t with a consequent increase in the 
torsional natural frequency p 2 . However, in behavior which is very 
similar to that noted in the preceding section fox- kt « 0 , increasing 
fi. causes an even more marked increase in k^., p l# and p t Thio 
result can be explained by noting that k t is ah (approximately) 1 Mv.s&x 
function of the tensile force .in the cab 34 ? whioh in propor- 

tional to T7 2 . These comparatively large values of kt also couple 
wobble and torsional motion. Evidence of this coupling is demon- 
strated in figure 5 for uncontrolled space station motion. 


o. 1/2 Deployed Configuration, fl « 4 and 6.23 rpm. 

The results for these two cases are presented in Table 4.4(c), 
and continue to demonstrate (a) the comparative insensitivity to 
the stage of deployment, and (b) the marked dependence of controller 
performance upon the nominal spin-velocity $7. The simulation results 
for the ft = 6.23 rpm case are presented in figures 7. a through 7.e. 

An inspection of these results confirms the effectiveness of the 
controller in providing initial wobble damping rates. The compara- 
tively minor effect of the lightly damped (s 3 = .1364) torsional 
oscillations on wobble motion is evident in Figures 7. a and 7.b. 

d .v 1/4 Deployed ft = 4 and 8.3 rpm . 

The results for these two cases are presented if* Table 4.4(d) , and 
reflect basically the same control cited above for the 1/2 deployed 
configuration. 
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TABLE 4.4(b): 3/4 Deployed Configuration 
Case Is ft » 4 rpm. = .418 rad, /sec., k^ = 1.726 x 10 s ft. lb. /rad. 

Input Data: p 1 /ft = .229 , p 2 /ft - 1.48 
Control Data: t; 1 - 1.0 , ? 2 = .1557 , ? 3 = . 2857 

kj « 710.9 lb. see. /rad. , k 2 - -1.694 x 10 4 lb. sec. /rad. 

k g = -2.070 lb. sec. /rad. , k 4 = 4.796 x 10 3 lb./rad. 

k 5 = -3.271 lb. sec. /ft., k 6 = -1.993 lb. /ft. , k ? - -.0547 lb. /ft. sec. 

Closed-Loop Pole Location: 

ft*, fty w 0 /w 

real -.1108 -.1108 -.1108 -.0288 

imag . 0 . 0 . 3716 . 7027 

Transient Performance: 

■ 13 ‘ 8 lb ’ ' Mmax " 15 ‘ 7 ft ‘ 


Case 2 : ft = 4.85 rpm = .507 rad. /sec. , k t ’= 2.544 x 10 5 ft . lb . /rad . 

Input Data: Pj/ft = .278 , p 2 /ft - 1.79 
Control Data: 1.0 , £ 2 = .1556 , £ 3 = .1667 

kj = 3079. lb. sec. /rad. , k 2 = -6.105 x 10 4 lb . sec. /rad. 
k 3 = -3.385 lb. sec. /rad. , k 4 = 3.750 x 10 3 lb./rad. 

k g = -1.247 lb. sec. /ft., k 6 = -8.198 lb. /ft., k ? = -.2367 lb. /ft. sec. 
Closed-Loop Pole Location 


^x ' fi y 


real 


-.1344 


xmag. 


0.0 


Transient Performance: 


u 


'max 


25.4 lb., w 


max 


w 0 / w 

-.1344 -.1344 -,0288 

.8530 .79 50 

= 16.0 f i . 
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TABLE 4.4(c): 1/2 Deployed Configuration 

Case 1: ST ~ 4 rpm. - .418 rad. /sec. , k t « 2.254 x xQ s f t. lb. /rad. 
l'nput Data: p : /U - . 229 , p 2 /H - 1.6 8 
Control Data: - 1.0, C 2 = .1364, r >3 ~ .3333 

kj = 373.8 lb. sec. /rad. , k 2 « -1.227 x 10 4 lb. sec. /rad. 
k 3 - 1.491 x 10 3 lb. sec. /rad. , k 4 = 4.691 x 10 3 lb./rad. 
k 5 = -3.578 lb. sec. /ft., k 6 = -1.621 lb./ft. ,k 7 = -.0402 lb. /ft. sec 
Closed-Loop Pole Location 

t w 0 /w 

real -.1107 -.1107 -.1107 -.0288 

imag • 0 . 0 . 3130 .8036 

Transient Performance: 

Mmax " 12a lbs - ' 1 w ' max = 14 ' 9 £t ‘ 


Case 2: ft = 6.23 rpm = .65 2 rad. /sec., k t = o.4 8 x 10 5 ft . lb . /rad. 

Input Data: p x /fi = .378 , P 2 /^ - 2.77 
•Control Data: 1.0 , £ 2 = ,1364, C 3 ~ *40 

k x = -174.96 lb. sec. /rad. , k 2 = -1.276 x 10 4 lb. sec. /rad. 
k = 6.907 x 10 3 lb.sec./rad. , k. = 1.181 x 10 4 lb./rad. 

3 ’ 

k s -5.928 lb. sec /ft., k 6 = -2.893 lb./ft. ,k 7 =-. 0685 lb. /ft. sec. 
Closed-Loop Pole Location 

P x , fiy w 0 / w 

real -.1726 -.1726 -.1726 -.0288 
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TABLE 4.4(d): 1/4 Deployed Configuration 

Case 1: n « 4 rpm. = .418 rad. /sec. , k t = 2.836 x 10 ? £t, Ib,/tad. 

Input Data: Pj/ft « .228 , p 2 /jT 83 2.52 
Control Data: 1.0, C 2 ~ . 1046 , ? 3 - .20 

kj « 2.156 x 10 3 lb. sec. /rad. , k 2 = “3.437 x 10 4 lb.sec./rad. 
k 3 « “1.2608 x 10 4 lb. sec. /rad. , k 4 « 3.372 x 10 3 lb,/rad. 
k 5 « -1.255 lb. sec. /ft., k s = -3.909 lb. /ft., k 7 = -.1113 lb. /ft. sec. 
Closed-Loop Pole Location: 

ft x , fiy w 0 /w 

real -.1104 -.1104 -.1104 -.02878 

imag . 0 . 0 . 5411 1.050 _ 

Transient Performance 

|ul v * 15.8 lb., |wl v = 16.5 ft. 

' 1 max ' 1 ' max 


Case 2 : S2 - 8.3 rpm. = . 867 rad. /sec. , k.^ - 1. 550 x IQ 6 ft. lb. /rad . 

Input Data: p 1 /'^ ~ .254 , p 2 /?J = 2.52 

Control Data: % l - 1.0, t, z ~ . 1046 , C 3 = 0.50 

kj » -712.6 lb. sec. /rad. , k 2 = -1.010 x x0 4 lb. sec. /rad. 
k 3 « 1.782 x 10 4 lb. sec ./rad. , k 4 - 1.892 x 10 4 lb. /rad. 
k 5 « -8.685 lb. sec, /ft., k 6 = -4.334 lb. /ft. , k 7 = -.0766 lb. /ft. sec. 
Closed-Loop Pole Location 

w 8 /w 

real -.2291 -.2291 -.2291 -.0288 
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CHAPTER V 

SUMMARY, CONCLUSIONS, AND EXTENSIONS 
5.1 SUMMARY AND CONCLUSIONS 

A dynamic model for cable-connected artificial-g space stations 
with an attached mcvable-mass-controller (MMC) system has been 
derived, and a linearized model for this system has been developed 
and its validity confirmed. This model has been used to examine 
the effectiveness of the MMC system in providing attitude stabili- 
zation for a space-station configuration reported by McDonnell- 
Douglas in Reference 8. The influence of the following three 
specific parameters were investigated: 

(a) Stage of deployment. Undeployed, 1/4, 1/2, 3/4, and fully 
deployed stages were examined. 

(b) Nominal spin-velocity U. This parameter was varied between 

4 rpm and an upper value depending upon the stage of deployment. 
(c> Torsional stiffness k t . This parameter was varied from zero 

(corresponding to a single cable) to an upper value calculated 
for the McDonnell-Douglas design. 

The MMC was required to provide critical wobble damping rates 
for all configurations without violating the peak-force-magnitude 
constraint, |u| 5 ^5 lb., and the peak-deflection constraint, 

| w | - 16.5 ft. As discussed in Section 3.6, damping for the control 
mass and the torsional motion was provided to yield a uniform 
damping constant for wobble motion, control mass motion, and rela- 
tive torsional motion, i.e., 

2C.P, = 2; 2 p 2 = 2c 3 p s 

where the are modal damping factors, and p x , p 2 , p 3 are the 




"natural frequencies" assoc: atod with wobble, torsional, and con- 
trol mass motion, respectiwly, In all cases £j was required to 
be unity (critical wobble damping)? hence, 

/ P 5 (37) 

The natural frequency, p 3 wa? positioned arbitrarily (within 
reason) for each case by a proper selection of the feedback gains? 
however, the wobble frequency p 2 and torsional natural frequency 
p 2 are physical parameters which depend on the configuration, (i.e,, 
spin velocity, asymmetry, torsional stiffness, etc.)* 

The procedure for comparing the controller's performance 
for the various cases is markedly simplified by the synthesis pro- 
cedure employed, since the "quality" of control is the same for 
all cases. Specifically, for all cases, wobble motion is critical- 
ly damped (Cj^l) ," and ‘all 'modal damping constants are equal. Hence, 
the significant point in a comparison of cases is the required 

maximum controi force |u| and control mass deflection I w I 

As previously noted, these two variables are of interest be- 
cause of the physical constraint on |w| , and the prescribed con- 
straint on |u| . They also indicate whether control torques arise 
primarily due to the control force u or control mass displacements 
w. Control torques directly caused by u are generated via the 
moment terms ur x and ury, while torques due to w are generated via 
product of inertia terms m c r x wfi 2 and m c rywTT 2 . From a controls 
mechanization viewpoint, the variable jul should be minimized, 

while I w I is maximized, i.e., Iwl - 16.5 ft. In other words, 

1 1 max 'ii max t 

one would like to employ 5 , a comparatively small and slowly-varying 
control force, and generate the required control torques primarily 


mi m ■ 


mm 


’■"iiTTria i iip— 1 


■ 
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through the "leisurely" manipulation of w. The feasi' ility of 
achieving this desired control character depends primarily upon 
the wobble frequency pj. As p* is increased/ the response of w 
to an error signal containing and tty is diminished/ and the 
required control torques can only bo supplied by increasing the 
magnitude of u. 

A similar (although lesser) problem arises when the torsional 
natural frequency p z is increased* However, the control logic 
which specifies ? 2 « P 1 /P 2 largely eliminiates any difficulty in 
controlling the torsional motion by simply reducing the required 
damping factor as p 2 is increased. The following points con- 
cerning this procedure should be noted: 

(a) As p 2 increases, the amplitudes of undamped torsional oscilla- 
tions are reduced, and the influence of this mode of oscilla- 

4 tion on the artificial-g environment in the crewquarters is 
reduced. Hence, the necessity for damping this mode of os- 
cillation is reduced as p 2 increased. 

(b) Attempts in this study to provide comp' rati vely large damping 
factors (i.e., ,4 f ? 2 1 .6) for configurations having com- 
paratively high values of p 2 were not notably successful in 
that a large control force (u 1 40 lbs.) was required, and it 
was generally necessary to reduce the wobble damping rates 
(.4 $ (, 5 .6). 

The results of all cases for I u I andjwl „„ are summarized 

in Table 5.1 and 5.2. Table 5.1 contain® the nominal (ST =4rpm.) 
cases, for both the zero and finite torsional stiffness cases. 

The entries of 0 and 1 for deployment stages indicate the non- 
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deployed and fully-deployed cases, respectively . The overspin 
cases arc summarized in TABLES 5.2. 

TABLE 5.1(a): Transient performance for Ti *= 4 rpm., and kfc * 0, 


Deployment 

Stage 

1 blmax 
(lb.) 

'tkf? 

0 

15.9 

1C. 2 

1/4 

19.0 

16.5 

1/2 

9.7 

16.1 

3/4 

19.1 

16.5 

1 

14.2 

16.5 


TABLE 5.1(b): Transient performance 

for sT « 4 rpm., and k t > 0. 

Deployment 

Stage 

'to 

lwlmax 
(ft. ) 

0 

15.9 

16.2 

1/4 

15.8 

16.5 

1/2 

12.1 

14.9 

3/4 

13.8 

15.7 

1 

.16.1 

16.3 


TABLE 5.2(a): Transient performance for overspin cases with k-^ - 0 . 


Deployment 

Stage 

S2 (rpm. ) 

1 u 1 max 
(lb.) 


0 

11.5 

24.4 

10.8 

1/4 

8.3 

19.0 

16.5 

1/2 

6.23 

34.0 

16.5 

3/4 

4.85 

28.0 

16.5 

1 

4.0 

14.2 

16.5 


TABLE 5.2(b): Transient Performance for overspin cases with k t > 0. 


Deployment ft (rpm.) ! u lmax l w lmax 

Stage • (lb. ) (ft. ) 


0 

11.5 

24.4 

10.8 

1/4 

8.3 

34.9 

9.7 

1/2 

6.23 

27.2 

13.7 

3/4 

4.85 

25.4 

16.0 

1 

4.0 

16.1 

16.3 

t 
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The results of these tables confirm the general conclusion 
of this study, namely, the MMC system can provide effective and 
generally satisfactory attitude stabilization for a wide range of 
space station configurations. In addition, the conclusions con- 
cerning specific parameter variations are as follows : 

(a) Stage of Deployment. For a fixod spin velocity sT, the changes 
in the' stage of deployment have a comparatively minor direct in- 
influence on controller performance. This statement can be con- 
firmed by inspecting Tables 5.1(a) and (b) . 

(b) Nominal spin velocity sT. This is the single most important 
factor on system dynamics. In general, effective control is much 
more easily achieved for the nominal spin-velocity H « 4 rpm. 
than the higher spin rates, which would suggest the advisability 
of spin-up or spin-down firings at the intermediate stages of 
deployment. 

(c) Torsional stiffness kf In general, the MMC was equally 
effective for the k t = 0 and McDonnell-Douglas designs. When k t 
is zero, the torsional motion is uncontrollable, but has no effect 
on the crew quarter artificial-g environment. Furthermore, the 
gains k t computed for kt - 0 proved to be generally satisfactory 
for k t up to 2000 ft. lb. /rad. For the finite torsional stiffness 
cases, 0 motion was coupled into the crew quarter"*, but the finite 
values of kt enabled effective control of 0. During the course of 
this study, preliminary values for kt were examined which were 
(erroneously) an order of magnitude larger than those cited in 
Section 4.4. However, gains were computed for these values (at 
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‘si « 4 rpm for tho fully-deployed configuration) which yielded 
an extremely effective control logic. At nose elevated values 
the torsional frequency p 2 was much larger than the wobble fre- 
quency p x , which in turn decreased and approached the rigid-body 
value (ab)'*’ / ^. As noted previously, the control logic used in 
this study minimizes any control difficulties associated with 
large values for p 2 . 

In some cases where comparatively high wobble frequencies pi 
were encountered, a smaller value of control mass rn c = 7 slugs 
was employed in an attempt to lower the peak force magnitude 
l u lmax* ^ e results were generally unsatisfactory, since this 
led to difficulties in satisfying the deflection constraint 
j w | <16.5 ft. 

5.2 EXTENSIONS 

A logical sequel to this investigation would consist of 
a similar study with an improved dynamics modeu. (such as that pro- 
sed by Nixon ) for cable-connected configurations. An improved ly- 
namics model would include the complete rotational degrees of free- 
dom for crewquarters and' counterweight, but would not’ require .extensiona- 
motion between the bodies . The. control logic and synthesis pro- 
cedure employed in this study should apply for an improved model, 
since basically only an increase in dimensionality is involved. 

Several practical problems which would arise in the imple- 
mentation of an MMC system remain unresolved. For example, the 
proper technique for measuring system variaoles (0, 6, etc.) is 
not clear. Furthermore, while the physical requirements of an MMC 
system would seem to be modest, a hardware investigation with the i 

f L: 


G4 


objective of fabricating and testing an MMC tc establish its 
operating characteristics would be of considerable value. 
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APPENDIX A: Closed-Xjoop Chftx avtcristic Polynomial 

The closed-loop polynomial, which is obtained by an ex- 

. i t i 

pansxon of the determinant is l~A~bk | , is of the form 

s* + a 6 s 6 + a g s s + a^s 1 ' + a 3 s 3 + a^s 2 + ajS l +a 0 = 0 (A.l) 

where each of the a.'s is a linear function of the unknown 
* 1 

gains i.e. , 

6 

a. = c-j + I a. .kn (A. 2) 

1 * - XJ J 

3=1 J 


From Eg. 

( 26 ) , 

the a 

. . coefficients 
13 

were determined as 

follows : 

a 6 : 

ct 6l 

■ w- 

a 

6 3 ~ 6 

2 

<V I V 




a. , 

6 4 

= a 6 6 

= 0, 

a 

65 = -d 

+ e 

: )/™c, a 

6 7 ~ ^ 


a 5 T 

“si 

= ft a . 
6 

3 '[a 

+ 

e (1 - b 2 

)] 




‘ 

a 5 2 

= ft 

0 

i b l 

t 

a s3 = a 

6> (b 

2” b i> 

• 



U 5 4 

= “6 3 

' a 5 

5 

“ 0 ^ s6 

= C V 

5 f0t 57 = 

0 


0i 4 i 

a 4 1 | 

= “SI 

[^ 2 (1 

+ 

b ? _) + P 






a 42 

= a 03 

[-ft 2 

(1 

+ b 2 + 

eb 2 

- eb 2 ) - 

(k. /I 2 
t' y 

)] 


a 4 3 

= ft 2 a 

6 3 ^ 

t 

ab 2 ) 






a 4 4 

= n a 6 

1 ^ b 2 

- 

b i> 






<*4 5 

= - (l/m c + r 

a 

y 

6 1 >Pi 2 







+ (ft 2 

/”c H 

e (1-b 2 ) (b 2 

- bi 

) - abj 

- b 2 ) ) 



“ r y a 6 1 + h Z ) 

- ft 2 r x 0^3 {a + e (1- b z ) z ~ b 2 - k t /(I 2 ft 2 )} 
a 4 G = 0 t <*47 = a 56 


M 
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a 


a. 


«i 


a, 


«3l 

1 B 

^ a 5 3 

2 [a (l+b2 } + e (l-b2) ] + 

(ak^/I 

2 )} 

y 

a 3 2 

ss 

0 ( g j 

{ 5 2 b x ( l+b 2 ) 

+ k t (b z /I 

1 

y 

+ b i/I 

y )} 

a 3 3 

= 


' a 34 = 

01 4 3 

/ a 3 s ~ 

0 



a 36 


a 4 s + 

2 (l/m c + 

r y a 

6i>V 


* 


a 3 7 

= 

«46 = 

0 , 






a 2 1 

= 

n 2 a 61 

(ft 2 b 2 + 

P 2 
1 1 

) 




a 2 2 

= 

^ a 6 3 

(ft 2 [eb 2 ( 1 “ 

b 2 ) ~ b 2 1 

- 

(.’ t/1 2 

y)J 

<*2 3 

= 

n 4 a 63 

ab 2 / 

a 2 4 

= a 3 3 ' 

a. 

!6 - 0 


^2 S 

= 

fi 2 a ( 1 - 

b 2 ) r x 

tt G 3 






* 

-O 2 (a/m 

c + r y a e 

a ) {k t (b 2 /I 1 y 

+ 

bj/I a y 

)+ S*bib*} 

a 2 7 

= 

a 3 6 







«i j 

- 

ft 3 aa 6 

3 (^ 2 b 

2 + 





a i 2 

= 

fi 3 a 6 x 

{ fi 2 b 1 b 

2 + 

^ t ^ b 2/ i * y 

+ 

b 1 /X Z y 

)} 

a i 3 

= 

0 , a 

1 4 = 

4 ab 2 

a 63 ' ^15 

: 

= 0 

* 

a i G 

as 

a 2S ' 

a !7 = ° 






“01 

= 

a 0 2 = 

^0 3 = 01 

0 4 5=1 

a 0 5 = a 06 

= 

= 0, a 0 

7 = a . 16 


The cj_ coefficients are: 


c 6 = c 4 = c 2 = 0 


c 5 = -fi 2 {e (l-b 2 ) (b 2 - b 2 ) + e/} 


cj = fl 2 ab 2 (k t /I ! y)+ (fi 2 b 2 + k t /I 2 y ) (fi 2 ab 1 - e^ 
-m^MbjryOtg i + ar x a ( 

-m c fi 2 r y a 6l k t /I l y 


6 3 


r X e < 1 - b 2 ) 2a 6 3 } 


c i ~ m c ^ H Ub 2 » 2 + k t /I 2 ) (b 1 r y a 6l + ar x a g3 ) 
+r y b 2 a 6 1 k t /I y } 




